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Abstract 

The Self Consistent Random Phase Approximation (SCRPA) and a Direct Analytical (DA) 
method are proposed to solve the Extended Hubbard Model in ID. We have considered an Extended 
Hubbard Model (EHM) including on-site and off-site interactions for closed chains in one dimension 
with periodic boundary conditions. The comparison of the SCRPA results with ones obtained by a 
Direct Analytical approach shows that the SCRPA treats the problem of these closed chains with 
a rigorous manner. The analysis of the nearest-neighbour repulsion effect on the dynamics of our 
closed chains shows that this repulsive interaction between the electrons of the neighbouring atoms 
induces supplementary conductivity, since, the SCRPA energy gap vanishes when these closed chains 
are governed by a strong repulsive on-site interaction and intermediate nearest-neighbour repulsion. 

Keywords: Extended Hubbard Model, SCRPA, energy gap 

Pacs numbers: 71.10, -W.75.10.Jm, 72.15.Nj 

1 Introduction 

The Hubbard Model is one of the most important models in the study of strongly correlated systems 
[1, 2]. In spite of its simple definition, it's believed to exhibit various interesting phenomena including 
metal-insulator transition [3, 4, 5, 6, 7], antiferromagnetism [8, 9, 10, 11], and superconductivity [12, 13]. 
In its description of interacting electrons on a lattice, the Usual Hubbard Model exhibits the competition 
between the usual kinetic and the on-site Coulomb interaction. In order, to establish a solution of the 
self-consistent equations taking into account the correlations of type 2p-2h (two particles- two holes), 
Jemai et al [14] proposed the Self Consistent Random Phase Approximation (SCRPA). They used RPA 
method [15, 16], and introduced the RPA excitation operators for pair: particlc-holc (p-h), where the 
vacuum represents the ground state of the system. These excitation operators are constructed from 
the linear combinations of creation and annihilation operators of pair: p-h. The quality of the SCRPA 
method to solve the Usual Hubbard Model has been investigated in a previous work by Ref. [14], in 
which, the authors showed the remarkable agreement between the SCRPA and exact results. 
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To explain other physical phenomena observed in different areas of the solid state physics like mag- 
netic and transport properties, an Extended Hubbard Model is proposed [17, 18, 19] taking into account 
the off-site interaction, with supplementary parameter: V (off-site interaction energy). 

In this paper, wc have analyzed the off-site interaction effect on the some local properties of the 
ID closed chains as: ground state energy, occupation numbers and energy gap. The paper is organized 
as follows. In sec. 2, we present the theory and describe, briefly, the SCRPA technique. In sec. 3, 
we present the SCRPA results and compare them with ones obtained by a Direct Analytical approach. 
Finally, we discuss the V effect on the dynamics of system. 

2 Theory 
2.1 Model 

The Extended Hubbard Hamiltonian is given by [17]: 

The first term of the Eqs. (1) represents the kinetic energy of electrons, where each electron has a 
possibility of hopping between different lattice sites. Cj ,j is the annihilation operator of the electron at 
a lattice site j with spin index a. g. is the creation operator of the electron at a lattice site i, so 
is the hopping integral from the site j to the site i. The second term represents the on-site coulomb 
interaction with energy [/, where rij^o- is the number operator of electrons at the site i with spin a 
. describes the effective off-site coloumb interaction between the electrons in the lattice sites i 

and j, with spin a and (j'respectively. The model (1) cannot be solved in a general case. In this work 
we will limit ourselves to a simplest case. We consider closed chains in one dimension with periodic 
boundary conditions, where each closed chain is organized alternately of two types of atoms noted 
and "_B" occupying the site i and i + 1, respectively, and characterised by three parameters: t (the first 
nearest- neighbour hopping), U (the on-site interaction energy) and V (the off-site interaction energy). 
With U > V > since the coulomb interaction is always repulsive and decreases with the distance. 

In the following of this paper, we suppose that U >- 2V. In this case, the ID highly correlated metal 
have an anti- ferromagnetic order [20]. Thus, we have (nj,crni±i,o-) — 0. Under these conditions, the 
Hamiltonian of our physical system is: 

Hii = -t ^(c-^^Cj+i,^ + c\j^i^aCi,a) + U (ni,|ni,| + ni+i,|ni+i,|) + ^ ^ n^^^nj+i (2) 

(7 (7 

In the first Brillouin zone — tt < A; -< tt, with the boundary conditions, the possible wave numbers 
are: fci = and ^2 = — tt. Thus, the Hamiltonian of our system can be developed in Hartree-Fock {HF) 
approximation as: 

Hhf = Ehf + XI {^1 {} ~ "■ki,aO-ki,a) + £2al^^„ak^,a^ (3) 

(T 

^ and ak,a are, respectively, the creation and annihilation operators of electrons with momenta k 
and spin a. 

The Eqs. (3) shows that, in the HF approximation, the system has two possible excitation energies: 
€i = t— ^^^Y^ ; €2 =t + ground state energy is 



Ehf = {HF\ Hhf \HF) = -2t+ -^^-^ 

Where \HF) = a\^^^a\^ \vacuum) is the HF ground state of system. Whereas the HF first exited 
state is taken as: \HF)* = aj.^i.a|.^ ^ \vacuum). Thus, we can define the HF quasiparticle operators by: 
&i,cr = Ofei cr ^"^^ ^2,0- = o-k2,<T- Therefore, we have hk,cr \HF) = for all k. 

Finally , using the usual Fourier transformation for the operator Ci^„, we can developed the Hamil- 
tonian Hn as function of the operators 61, o- and 62,(7- 

H = Hhf + Hk=o + Hk=-^ (4) 

where 

U + V 

Hk^o = — - — (nfc2,t - rifcjj) {nk2,i " "fci.i) 

Hk=-^ = (jf + + jI) 

with 

Hk=o and Hk^_T^ take into account, respectively, the correlation between the number operators of 
the type: nki^afi'kj,a' and between the magnetic moment operators of the type: JI-J~, ■ 

2.2 SCRPA approach 

The Random Phase Approximation (RPA) [15], was used to solve the Usual Hubbard Model [14, 21]. 
The RPA is an approach which treats seriously the correlations of system, and attempt to minimise the 
system energy. In order to apply the formalism of SCRPA to the Hubbard Model, it is convenient to 
use the ph-RPA regrouping, which regroup the physical system on pairs: particle- hole (p-h). We can 
then, define the RPA excitation operator as: 

QI = JZ'^^lhblbi - y^hbhbp) (5) 

Where h (and p) are the momentum below (and above) the Fermi level. The Eqs. (5) shows that the 
excitation in the ph-RPA is done only by the creation or (annihilation) of pair: particle-hole via the 
operator bpb^f^ ipkbp) with the amplitude a;^^ iVph)- "^^^ corresponding excited state of this excitation 
operator is \v) = \RPA), and the corresponding excitation energy is: 



{RPA\ Q„ H,qI \RPA) 



{RPA\ [Q„,Oy \RPA) 
Where \RPA) is the vacuum of this RPA excitation operator: 



(6) 



g„ \RPA) = (7) 

We consider that this vacuum represents the ground state of the system. It's clear, that the RPA 
approach will be more rigorous than the HF approximation in the treatment of the correlation of 
system. Since, its ground state is not defined explicitly like the HF ground state, but, it is defined only 
via the condition of the Eqs. (7). Thus, the RPA ground state can contain the correlations. 



The minimization of leads to usual RPA equations of type: 



/ A B \ I x" \ _ I \ 
V -B* -A* )\y- )^ [y" J 
Where A and B are two square sub-matrix. 

For the charge and spin longitudinal sector (5 = 1, mg = 0) [14] , we consider only the RPA excitation 
operators which conserve the spin. Where the excitation is done only by the creation or annihilation of 
the pair: particle-hole with same spin. 



(8) 



with 



Ma = ni,<7 + n2,a 



In this sector, the sub-matrix A and B take the form: 



A = 



Hi 



Where the RPA matrix elements are given by: 



1 



and 



and 



V(l-M<,)(l-M,,) 



B = 



Hi I, J a 



B'^'^ B 
-Bit B 



{[j-,[Hn,J-]]) 



With the relations of the orthonormality conditions of the set |(3v;<3l'|, we have expressed the 
elements of A and B by the RPA-amplitudes, and therefore we have constructed a closed system of 
non linear coupled equations, which we have solved numerically by iteration leading to the SCRPA 
solutions. For more details, sec Rcf. [22]. In the present paper, we will compare these SCRPA solutions 
with ones obtained by a Direct Analytical approach for two local properties of our closed chain; the 

S C RP A 

SCRPA ground state energy: E^^ = {RPA\ Hu \RPA) and the occupation number of states with 
spin a and momenta k below the Fermi level : n^^^^ ^ = 1 — {RPA\ J^J^ {RPA). 

3 Results and discussion 

In this section, we present our SCRPA results and compare them with ones obtained by a Direct 
Analytical (DA) method. 



3.1 DA method 

To take into account the correlation of the first order for the ground state, we suppose that the Direct 
Analytical ground state wave function contains the correlation of type {2p — 2h) , therefore the ground 
state is: 

io) = (cs+cf^.j;j;)[7?F) 

The Hamiltonian of system in the subspace { \HF) ; j| ■ \IIF) } is: 



H= -''^"^ -^/M (9) 
\^ -/3/2 2t + a J 

Where 

a^U + V and f3^U -V 
This Hamiltonian has two roots: 

Eo = a 

and 

Ei = a-\ 

We take the eigenvector corresponding to Eq as the exact ground state, we have then: 

|0) = (cos(0) + sin(</)) j|' \HF) (10) 

with 

^ = arctan(- ^ ^ 



4t+V 0^ + 16^2 

Finally, the Direct Analytical gronnd state energy and the occupation number of states with spin a 
and momenta k below the Fermi level are calculated, respectively, via: 

E^J = (0| H |0) and = 1 - (0| J^J" |0) 

Using the Eqs. (9) and (10), we can find: 

Ks = -2i cos(20) + ^(1 - sin(2<^)) + ^(1 + sin(2<^)) (11) 

,a = cos^icj)) (12) 

Where, the parameter <p is given by: 

u-v 

(b = arctani . ) 

4i+ V(f/-F)2 + 16t2^ 

3.2 SCRPA and DA results 

In Fig.l (Fig.2), we have plotted the variation of E^^ and nk-ikF,a of SCRPA and DA methods as 
function of U/t for different values ofV/t . 

The Fig.l and Fig.2 show that the agreement between the SCRPA and DA methods is very good 
for any value of U/t and V/t. It is clear that the SCRPA method loses partially its precision for the 
high values of U/t, but as same on an important interval [0 ; 6t\ the method keeps all its precision. 
We can conclude, thus, that the SCRPA method is efficient to treat seriously the off-site interaction of 
type rii o-f^i+i.-cr- The curves show, also, that for U < V . we have no SCRPA solution. The program 
associated to SCRPA method begins to turn and give the solutions for U ^V; whereas the DA method 
gives solutions even for U ~< V. Thus, we can conclude that the SCRPA method is more realist than 
the DA method; knowing that for a real physical system, the on-site interaction is always stronger than 
the off-site interaction {U >- V). 



Figure 1: Ground state energy of SCRPA and DA methods as function of on-site interaction energy U jt 
for different values oiVjt. 



Figure 2: Occupation number of states with spin a and momenta k below the Fermi level of SCRPA 
and DA methods as function of on-site interaction energy JJ jt for different values of Vjt. 



Figure 3: Ground state energy of SCRPA and DA methods as function of off-site interaction energy 
V=t at U = 2t. 

In Fig. 3, we have plotted the variation of E^^ of SCRPA and DA methods as function of V/t at 
U = 2t. 

The Fig. 3 shows that the SCRPA and DA results are practically equal for any value of V/t. Noting 
that the program associated to the SCRPA stops to give the solutions, when V becomes the order of U 
(y ~ J7). This figure shows, also, that the ground state energy of our closed chains has, practically, a 
linear V dependence. 

In summary, we have compared the SCRPA results with ones obtained by a DA method, wc have 
found that SCRPA method solves the Extended Hubbard Model for closed chains in ID exactly for any 
value of U and V. 

3.3 Effect of V 

In order to analyse the effect of V on the dynamics of our closed chains, we define the energy gap Ae 
as the difference between the first excited state energy ei [22] and the ground state energy Escrpa- 

As = ei — Escrpa 

In Fig. 4, we plot the variation of the energy gap A^; as function of the repulsive on-site interaction 
energy U for different values of the off-site interaction energy V, with 2V -< U. 

First, we disregarded the off-site interaction: V = (Usual Hubbard Model case ), the corresponding 
curve shows that the gap energy A^; decreases when U increases. We deduce that the repulsive on-site 
interaction {U >- 0) increase the conductivity of system, since the repulsion between the two electrons 
of one site encourages every electron to jump to the neighbouring site. Then, we have introduced the 



Figure 4: Energy gap as function of on-site interaction energy U/t iov different values of V/t 



ofF-sitc interaction (V 7^ 0), and we have plotted the variation of the energy gap as function of the 
repulsive on-site interaction energy U for different values of the interaction off-site energy: V/t. The 
curves show that the effect of V on the energy gap is very weak when the on-site interaction is 
weak. For instance, at ~ 2f the effect of F on A^ is the order of 18% for an intermediate off-site 
interaction {V/t = 1). Whereas, when the on-site interaction energy increases, the effect of this off-site 
interaction will be more remarkable. For instance, at f/ ~ 6i (strong on-site interaction ) the effect of 
this intermediate off-site interaction {V/t = 1) on A^; will be the order of 50%. The curves show also 
that when U increases, the gap energy A^; decreases and tend to a constant value Ao at critical on-site 
interaction energy Uc- This limit value Ao depends on off-site interaction energy V, noticing that Aq 
vanishes for the intermediate repulsive off-site interaction {V/t = 1). We deduce that the repulsive off- 
site interaction increases the conductivity of system, especially, in the case of an intermediate repulsive 
off-site interaction, the gap between the first excited state and the ground state energy vanishes, and 
therefore each electron can jump between different sites without any loss of energy. 

4 Conclusion 

In this paper, we have used the SCRPA and DA methods to solve the ID Extended Hubbard Model in 
order to study some local properties of closed chains in one dimension with periodic boundary conditions 
as: ground state energy, energy gap and occupation numbers. First, We have compared the SCRPA 
results with ones obtained by the Direct Analytical method; we have shown that the SCRPA method 
solves the Extended Hubbard Model for finite number of sites exactly, for any value of U (on-site 
interaction energy) and V (off-site interaction energy). Then, we have analysed the V effect on the 
dynamics of system. We have shown that the energy gap decreases with V. Particularly, we have found 



that with an intermediate repulsive off-site interaction {V/t = 1) this gap between the first excited state 
energy and the ground state energy vanishes, and therefore the electrons can jump between different 
sites without any loss of energy. 
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